NIST-GCR-93-634 


FIRE INDUCED THERMAL FIELDS IN 
WINDOW GLASS I - THEORY 


TILE COPY 
A. A. Joshi and P. J. Pagni D0 NaT REMOVE 


University of California 
Berkeley, CA 94720 


NST 


United States Department of Commerce 
Technology Administration 
National Institute of Standards and Technology 


NIST-GCR-93-634 


FIRE INDUCED THERMAL FIELDS IN 
WINDOW GLASS I - THEORY 


A. A. Joshi and P. J. Pagni 
University of California 
Berkeley, CA 94720 


December 1993 


* vert OF 


Sponsored by: 

U.S. Department of Commerce 

Ronald H. Brown, Secretary 

Technology Administration 

Mary L. Good, Under Secretary for Technology 
National Institute of Standards and Technology 
Arati Prabhakar, Director 


Notice 


This report was prepared for the Building and Fire 
Research Laboratory of the National Institute of 
Standards and Technology under grant number 60NANB1D1168. 
The statements and conclusions contained in this report 
are those of the authors and do not necessarily reflect 
the views of the National Institute of Standards and 
Technology or the Building and Fire Research Laboratory. 


abe 


Fire Induced Thermal Fields in Window Glass I - Theory 


A. A. Joshi* and P. J. Pagni 
Mechanical Engineering Department 
University of California at Berkeley 
Berkeley, CA 94720 


Abstract 


Window glass breaking plays an important role in compartment fire dynamics as the win- 
dow acts as a wall before breaking and as a vent after breaking. Previous work suggested a 
model for the time to breakage of a window glass exposed to a particular fire. In this paper, the 
glass thermal fields obtained using that model are examined in detail. The temperature field 
dependence on heat transfer coefficients, radiative decay length and flame radiation is explored. 
The results show that the glass surface temperature increases with a decrease in the decay length 
and increases with an increase in flame radiation heat flux. Early in the fire, the glass temperature 
may be higher than the hot layer temperature due to direct impingement of flame radiation. Later 
the glass temperature lags the hot layer temperature. The variation of the time to breakage as a 
function of the shading width and decay length is also presented and the results indicate that the 
ereaking time decreases with an increase in the shading width and decreases with a decrease in 
, decay length. Heat flux maps for typical conditions indicate that most of the heat influx is stored 


in the glass, increasing its temperature. 


* Present Address: Mechanical Engineering Department, University of California, Santa Barbara, CA 93106; 
to whom correspondence should be addressed 


Nomenclature 


A-G __ constants 


Bi Biot number, AL /k 


Cc specific heat capacity 

E Young’s modulus 

F kemels 

H half-width of the window 
h heat transfer coefficient 

I radiant heat flux directly from the flame 
j(t)  =I(t)L/kT,, see Eq. (6) 
k thermal conductivity 

L glass thickness 

l decay length 

Q energy 

q dimensional heat flux 

S shaded length 

T temperature 

u dummy variable 
subscripts 


1 ambient (outside) side of glass pane 
2 compartment side of glass pane 

b breakage 

c characteristic, convection 


f flame radiation 


initial 


r radiation 
‘Ss stored in glass 
y _ diffusion in shaded portion 
co ambient 
Greek Characters 
a thermal diffusivity 
6 thermal expansion coefficient 
Y dimensionless decay length, //L 
€ emissivity 
$ dimensionless heat flux 
6 dimensionless temperature 
t dimensionless time 
€& dimensionless x coordinate 
oO _ stress 


INTRODUCTION 


Window glass breaking in fires is an important practical problem, since a window acts as a 
. wall before breaking and as a vent after breaking. Thus the knowledge of the time to window 
breakage is crucial in order to predict the evolving compartment fire dynamics. Emmons [1] 
highlighted the importance of glass breaking in fires, and, based on Harvard experiments [2], sug- 
gested thermally induced tensile stresses as the mechanism for glass breakage. All window glass 
has its surrounding edge covered by an opaque frame or gasket. Glass is a relatively poor thermal 
conductor. The central portion of the window is heated by infrared radiation and hot gas convec- 
tion, while the frame protects the shaded edge so that it remains near its initial temperature. The 
thermal expansion of the central portion places the covered edge in tension until it cracks. These 
cracks then bifurcate and quickly propagate causing the window to become a vent. This scenario 
suggests a simple criterion [3] for the temperature rise required to break glass windows in fires 


which follows directly from Hooke’s law [4]. 
BAT =0,/E : (1) 


where AT is the temperature difference between the glass edge and center, B is the thermal 
coefficient of linear expansion, GO» is the breaking stress and E is the Young’s modulus for glass. 
Keski-Rahkonen did excellent analyses [5,6] of the limiting cases where: 1.) there are no tem- 
perature gradients across the glass thickness so a constant heat input can be treated as a 
volumetric heat source; 2.) convective and linearized radiative heat losses are constant and identi- 
cal on both sides of the glass so it can also be treated as a volumetric heat sink and 3.) all gas 
temperatures, the initial temperature and the temperature of the outer edge are constant and ident- 
ical. His analytic results confirm Eq. (1) as a lower limit for AT. While the restrictions make it 
difficult to apply these analyses directly to compartment fires, they do describe small-scale exper- 


iments with large, constant radiative heat influx [2]. 


Skelly, Roby and Beyler [7] did an experimental study of glass breaking in compartment 


fires. They burned liquid fuel in pans inside a compartment which had a glass window and 


recorded the temperature histories of the glass and hot gas layer inside the compartment. This 
experimental study also indicates that Eq. (1) is a reasonable first approximation for the glass 
temperature at breakage. What remains is to model energy transfer between the fire, the window 
glass and the ambience, to predict the time at which the glass breakage temperature is reached. 
Techniques to accurately calculate the history of the central glass temperature profile, T (x,t), by 
taking the Laplace transform of the temperature and numerically integrating the resulting equa- 
tions, have been described [8]. Two-dimensional temperature histories, T(x,y,t), and mean 


stress histories, 6,, (y ,t), were also presented. 


This paper examines in detail the thermal fields obtained from that model [8]. The follow- 
ing problems are addressed: 1.) the nature of kemels arising from Laplace Transform of the tem- 
perature; 2.) parameter effects on the temperature profiles within the glass thickness; 3.) the 
dependence of the temperature histories of the central exposed glass surfaces on the Biot number, 
the dimensionless decay length, y, and the unsteady dimensionless incident radiative flux, /(t); 
and 4.) the dependence of the dimensionless breaking time, t,, on the dimensionless shading 
width, s/L. The dependence on the heat transfer coefficient at the outer surface is not explored 


because it is not expected to vary significantly in a fire environment. 


THEORY 


Consider the window shown in Fig. 1. For large shading, s/L22, and fast heating, 
Cty /s2s1, where s is the width of the shading and L is the glass thickness, the temperature differ- 
ence AT in Eq. (1), corresponds to the temperature difference between the initial temperature of 
the unheated outer edge and the transient temperature of the uniformly heated central section of 
the glass [8]. The goal then is to determine the temperature in the large central section of the 
glass as a function of depth into the glass, x and time, t. Significant temperature gradients, 
OT /dx were shown to exist since the heat source is on the inside of the window and the sink is on 


the outside. The unshaded glass is uniformly heated, so d/dy and d/dz are zero. The goveming 


equation is 
2 -xil 
pc St =p ST 1), (2) 


where /(t¢) is the incident radiative flux directly from the fire which is at sufficiently short 
wavelengths that its distributed intemal absorption needs to be included [9] and / is the decay 


length in the glass. With the assumption that the glass is grey to other radiation, the initial and 


boundary conditions are 
att=0,T =7;, (3) 
etx =9, Kor =h2(T 200(t -T (0,t)) + €2e0OT 200°(t) — €0T 4(0,t) = g2(t), (4) 
atx=L, 49 =h,(T(L,t)-T.(t)) + €0T 4(L ,t) — €).0T j0*(t) = 9 1(t), (S) 


where side 1 is toward the outside ambience and side 2 is toward the hot layer inside the compart- 


ment. 


With the definitions 
E=2;t= Of! 9— T-Tj_. =9o /IEB-T;39 Be q\ >d2= q2 j@= I(t (6) 
Lite pr Pps Os pt Te = OH/EB— Tes O1= pri = pro I= ere 


the dimensionless governing equation is 


00 _ 0°70. «vet 
tan ee + j(T) y’ (7) 
with dimensionless initial and boundary conditions 
at t=0,6=0, ; (8) 
at &=0, -# = $,(1), (9) 
at &= 1- $= 100). (10) 


Explicit expressions for the dimensionless heat fluxes, 6;(t) and 2(t), are given in Appendix I. 


This equation is solved using a Laplace transform on time, 


9° = Joe? tar. 
The solution for both surface temperatures is given by 


a % G 
0(0,t) =2 J F ,(0,u),(t-u2)du +2 i F2(0,u)2(t-u2)du + 2 i F3(0,u)j (t—u2)du 


Y 


and 


7 ier 
6(1,t) =2 [Fidunene-uau +2 frau yenceuydu +2 [Fs(u)y(-u?)du, 


where the F; (E,u) are kernels which are functions of y, & and u and are 


for short time, u <0.4, 
FiGu)= + pT + die oe 


= ~G=2kP 
FE) = qr on on +h 4u2 ; 


and 


BPR ae 
: xe 408 Berd tas Ve ees cnr Meath V4}. 
+ Sere aust) 2] MPS ope Bester 5] pi 
. Le¥ Berd yet] OF are Bah) ae 
+ ord ENE) Ford eg] OF] 


(12) 


(13) 


(14) 


(15) 


(16) 


and for long time, u > 0.4, 


F,(€,u)=—u(1+2 pGily cos(k né)e~Hnu?), (17) 
F (Eu) u (142 XCD costk m(1—EyeP"*), (18) 
and 


1 
—aliey 2 S Cl) cosk —k2ny? 
F3(E,u)=ul e ¥}—y+ 7 yy Sate 2 


_| wa 2 < oes 1-€)) ,-k2ntu? 1 
ry A =k o ee 


The numerical procedure chosen was the trapezoidal rule with constant time steps and thus vari- 
able Au, since the 9; are functions of 8 which is known only at each time step. Newton- © 


Raphson’s method for finding roots of non-linear equations was also utilized. 


KERNELS 


Figure 2a shows the variation of F (0,4) and F (0,4) with u. From Eq. (12) it is evident 
that the area under the curve F ;(0,u) from u =0to Vt is equal to the temperature of the inner sur- 
face 6(0,t) for 6; =1/2 with 62=/ =0. Since 6; represents heat loss from the outer surface, 
6(0,t) should decrease with time when 9;>0. Thus F \(0,u) monotonically decreases with u. 
Also, since the initial condition is 80,0) =0, a positive >, implies that 6(0,t) should become 
negative. Therefore, F \(0,u) is negative. The magnitude of F;(0,u) remains close to 0 till u 
approximately reaches 0.4. The reason is that the thermal diffusion wave emerging at the outer 


surface due to the sudden heat loss at t=0 takes some time to reach the inner surface. 


The variation of F 2(0,u), also shown in Fig. 2a, is analogous to that of F ;(0,u) for ;=/=0 

and $2=1/2 except that now 2 represents heat gain at the inner surface when $2>0. Thus F 2(0,u ) 
_increases with u since the area under this curve from u=0 to Vt is equal to the temperature at the 
inner surface, 0(0,t). Also, the temperature of the inner surface instantaneously increases at t=0 


when 62>0. Therefore, F(0,u) begins with a positive value. This result can also be obtained by 
noting that the limit of F2(0,u) as u-0 is equal to fe: The plots of F;(1,u) and F(1,u) are 
not shown here since F ;(1,u) =—F2(0,u) and F2(1,u) =—F ,(0,u), 


Figure 2b shows the variation of = 3(0,u) with u for different values of the dimensionless 


decay length, y. As suggested by Eq. (12), the area under the curve, 4F 3(0,u) from u =0 to Vt 
is equal to 6(0,t), the dimensionless temperature of inner surface for j = 1/2 with >; = $2 =0. For 


a constant value of y, ie 3(0,u ) increases with u, since the imposed radiative flux increases the 


temperature of the inner surface, x=0. At a particular value of x, HF 3(0ut) decreases as Y 


increases because a smaller amount of net radiation is absorbed in the glass. In the limit as tT, 


the variation of F'3(0,u ) with respect to u is linear since the diffusion effects cease. 


The variation of 4F 3(1,u4) with u for different values of y is shown in Fig. 2c. The area 
under this curve from u =0 to Vt is equal to 6(1,u), the dimensionless temperature at the outer 
surface, when j = 1/2 and $;=$2=0. Here again, al 3(1,u) increases with u because the imposed 


radiative heat flux increases the temperature of the glass. In the limit as t--~, the variation of 
F3(1,u) with respect to u is linear due to on diffusion effects. Fig. 2c however, differs 
significantly from Fig. 2b in the variation with respect to y. This can be explained as follows. In 
Fig. 2c, two different time zones, u <0.4 and u>0.4 can be identified. The earlier time zone 
occurs when the thermal diffusion wave emerging at x=0, resulting from the temperature differ- 
ence between the surfaces, has not reached x=L. In this time zone, the temperature of the outer 


surface, Q(1,1), increases only due to the radiative flux. Therefore, as y increases, more radiation 


10 


reaches the outer surface and 6(1,7) and aF 3(0,u) increase. In the later time zone, the thermal 


wave emerging at x=0, reaches the outer surface x=L.. Here, 6(1,7) decreases with an increase in 


_y because a smaller amount of net radiation is absorbed by the glass. 


EXACT SOLUTION 


An exact solution to Eq. (7) subject to Eqs. (8-10) can be obtained for the case of linearized 
radiation and constant heat transfer coefficients, ambient temperatures and emissivities on both 
sides of glass. Here, the radiative heat loss from both surfaces of the glass is linearized giving 
effective heat transfer coefficients. Therefore all terms in Eq. (3-9) of the appendix which con- 
tain emissivities are dropped. The technique used to solve this problem is separation of variables. 
First, the solution is divided into a steady state and a transient part. The steady state part of the 
solution is then trivial and the transient part is obtained by solving the resulting eigenvalue prob- 


lem. The solution is 


6(€,t) =uE,t)+v), (20) 
where 
v(G)=p§+q, (21) 
and 
— — Bi1Bi2(02~ — 91.) _ Bi;01.0 + Bi2020 + Bi 1Bi202. 
Pe Bi, + Biz +Bi Biz and q = Bi, + Biz +Bi Biz . (22) 


The constants 02.., 0}... Bi; and Bi2 are defined as 


ge ide li” in oh di leorl Lte pe Bi ols 
02.0 = 1 » B10 = Bi; = —— and Bi2= ——-, (23) 


where h ; and h2 are the heat transfer coefficients. The solution for u(E,t) is given by 


‘ a2 Cae + fem, (n)dn| on), (24) 
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where the characteristic functions 6, are 
On (E) = An COS(An §) + Bi2sin(An §), (25) 
and the eigenvalues A, are obtained from the solutions to 


_ An? —BiBi2 


cota, = "h,, (BI Bip) (26) 
The constants C,, weight functions, w, and Pes N(A,) appearing in the equation for u (E,t) are 
defined as 

; =—yi Eon G48, wn (2) = sere (De, Od, @7) 
and 

NOnd= 4] On2 +01 | 1+ Pe ei. (28) 


Comparison between exact and numerical solutions showed that, for a dimensionless time 
step of 0.002, the agreement was within 0.5%. The exact solution however cannot be used to 


simulate real evolving fire environments. 


Profiles 


In this section, the temperature profiles obtained. in glass using the exact solution are 
explored as functions of different parameters to estimate the temperature distribution within the 
glass. Figure 3 shows the temperature profile as a function of the dimensionless hot layer tem- 
perature, 02... Here, 6; = 0, j =1, Bi; =0.1, Biz=0.5 and t=0.3. These values are chosen 
because they are typical of real fires. For example, for a glass window of 5 mm thickness, and 
thermal conductivity of 0.76 W/m-—K , Bi; =0.1 would correspond to Ay = 15.2 W/m?-K which 
is typical of natural convection in air on the outer surface exposed to ambience. The inner sur- 
face heat transfer coefficient would then be equal to 76 W/m2-K which is typical of forced con- 


vection by the hot fire gases and J = 28 kW/m2, which is typical of a compartment fire. For 


i be 


these parameters, the temperature of the glass increases with an increase in 62,. due to a larger 
heat influx. The maximum temperature of the glass occurs at the inner surface as long as the gas 
temperature remains higher than the glass temperature otherwise it occurs in the interior of the 


glass since the glass loses heat on both the sides due to convection. 


Figure 4 shows the temperature profiles as functions of j. Here 62,.=5, 0). = 0, T=0.3, 
Bi2=0.5 and Bi; =0.1. The temperature increases with an increase in j as expected. The varia- 
tion of temperature profiles with Biz showed that the temperature at the inner surface increases 
with an increase of Biz when 62..=5. This is a typical value of hot gas temperature in real fires 
at the time of breakage. In all the cases considered significant temperature gradients exist within 
the glass thickness implying that a uniform glass temperature approximation is not valid in real 


fire situations. 


PARAMETERS 


Figure 5 shows the temperature history of the inner glass surface, for different values of 
dimensionless decay length, y. Biot numbers are held constant and the hot gas temperature was 
calculated from FIRST [10] for a recent full scale fire [8]. Here J = 10 kW/m? and represents a 
window adjacent to the fire. The temperature of the inner surface increases with a decrease in Y 
because a small value of y implies that most of the incident radiation is absorbed within a short 
distance into the glass thickness thus increasing the temperature of the inner surface. For large 
values of y, less net radiant flux is absorbed in the glass and the glass behaves more like a tran- 
sparent medium. A similar calculation was carried out to observe the effect of Biot number, Biz 
on the inner surface temperature. Here / = 0 and the variation of the hot gas temperature was the 
same as in Fig. 5. The inner surface temperature increased with an increase in Biz due to a higher 


heat influx as expected. 


The temperature history of the inner glass surface with respect to varying input radiative 


flux function, j(t) was also explored holding all other parameters constant. The temperature 


1 


history of the inner surface was seen to roughly follow the integral with respect to time of the 


incoming radiation. 


TIME TO BREAKAGE 


The time to breakage is defined as the time at which the tensile stress at the shaded edge of 
the glass reaches the breaking stress. Fig. 6 shows the variation of the dimensionless time to 
breakage as a function of the dimensionless shaded width, s/L and the decay length, y with 01/62 
=0.33, 62 = 0.5 and j/o2 = 1.0. The thickness of the window glass is assumed to remain constant. 
The time to breakage increases with a decrease in shading thickness. The reason is that the ten- 
sile stress at the shaded edge is proportional to the difference between the exposed central and the 
shaded edge temperatures. If the shaded width is small, the edge temperature increases due to 
thermal diffusion, thus decreasing the tensile stress. Therefore, the central exposed surface tem- 
perature must be higher to reach breaking stress and hence the breaking time is longer. The © 
breaking time also increases with an increase in decay length because more radiation is transmit- 


ted through the glass as y increases. 


ENERGY AND HEAT FLUX MAPS 


In this section, the magnitudes of energy and heat flux from various inputs and outputs for 
the window glass are represented in terms of stacked bar charts. The heat inputs are convection 
and radiation from the gases, Q, and g, and Q, and q, respectively, and radiation from the fire, 
QO, and gf. The outputs are defined as heat oe from the outer surface, Q, and q}, radiative 
heat loss by emission from the inner surface, Q2, and q2,, heat diffusion into the shaded portion 
of the glass, Q, and q, and the rate of energy storage in glass, Q, and g,. The heat diffusion into 
' the shaded glass was estimated by assuming a hyperbolic tangent temperature profile for the tem- 


perature averaged from 0 to L under the frame. 
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Figure 7a shows the energy distribution evolution for the case of Fig. 5 where the glass is 
far away from the fire so that j(t)=0. The hot gas temperature variation was calculated using 
‘FIRST [10] for a recent full scale fire. Thus the glass heats up only by convection from hot gases. 
It is clear that very little heat diffuses into the shaded portion of the glass. This occurs, even 
though the heat flux is quite substantial (Fig. 7b) due to the very small cross sectional area of the 
window glass. The heat loss from the outer side (1) is also quite small as the temperature of the 
outer surface remains near the ambient temperature. The radiative heat loss from side 2 (the 
inner surface) is the largest loss since the glass temperature is a maximum at the inner surface. 
Most of the energy is stored in the glass increasing its temperature. The heat flux distribution 
evolution is similar to that of the energy distribution with the exception of the heat flow under the 
frame. The storage heat flux was calculated as the rate of change of internal energy integrated 


over the glass thickness. 


CONCLUSIONS 


The variation of the kernels obtained in the solution of the glass surface temperatures are 
explored as functions of time and decay length. The magnitudes of all kernels increase with time 


as the area under the kernels represents the temperature of the glass surfaces. The kemels 
4F 3(0,u) and 4F 3(1,4), which represent the temperatures of inner and outer surfaces respec- 


tively, when insulated to convection, decrease in magnitude as y increases because a smaller 
amount of radiation is absorbed in the glass. The dependence of temperature profiles on various 
parameters is also examined and it is concluded that in a fire environment, a uniform approxima- 
tion for the glass temperature distribution with depth is not valid due to the presence of large tem- 
perature differences between the inner and outer surfaces. The dependence of glass surface tem- 
perature histories on various parameters is investigated. The surface temperature increases with a 
decrease in the decay length, y. The temperature also increases with an increase in dimensionless 


heat transfer coefficient on the inner side, Biz, if the hot layer gas temperature is higher than the 


bes 


glass temperature. Also, the inner surface temperature varies roughly as the integral with respect 


to time of the unsteady radiative flux, j(t). 


The variation of the breaking time of glass windows is determined as a function of decay 
length, y, and window geometry, s/L and H/L. Breaking time decreases as decay length 
decreases, since more net radiation is absorbed. A large shading also decreases the breaking 
time. The magnitude of the energy and heat flux from different sources and sinks is compared in 
bar charts and it is observed that most of the heat is absorbed within the glass increasing its tem- 
perature. The most significant loss term is radiation from the interior glass surface. The heat 
diffusing into the shaded edge is quite small because the cross sectional area of the glass is small. 
The heat loss from the outer surface is low because the temperature difference between the outer 


glass surface and the ambience is usually low. 
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APPENDIX 
HEAT FLUX EXPRESSIONS 


The heat fluxes in Eqs. (9) and (10) are 

2(t) =A +B (0,1) + C 62(0,7) + D 63(0,1) + E 64(0,7) 
and 

$,(t) =F +G0(1,t) — C 6(1,1) — D 67(1,t) — E 64(1,1) 
where 


re hoL (T 2.0(t )-T;) fi E20OLT 0(t) _ eoLT;4 
Figuh kre. e, All DARRAKT 92 Bre 
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Figure Captions 


Figure 1. Window geometry showing the coordinate, x and length parameters. x is the depth into 


the window, s is the width of the shading, H is the half-height and L is the glass thickness. 


Figure 2(a). Variation of kernels F ,;(0,u) and F (0,4). The area under the curve F ,(0,u ) from 0 
to Vt is the temperature of the inner surface for the sample case, 9; = 1/2 and j =42=0. The area 
under the curve F2(0,u ) from 0 to Vt is the temperature of the inner surface for the sample case, 


$2 = 1/2 and j = 9; =0. 


Figure 2(b). Variation of a 3(0,u) with u and y. The area under this curve from 0 to Vt is the 


temperature of the inner surface for the sample case, j = 1/2 and 6; =$2=0. The symbols indi- 


cate the numerical data points at fixed u intervals (Au =0.01). 


Figure 2(c). Variation of yF (Lu) with u and y. The area under this curve from 0 to Vt is the 


temperature of the outer surface for the sample case, j = 1/2 and 6; =$2=0. The symbols indi- 


cate the numerical data points at fixed u intervals (Au =0.01). 


Figure 3. Temperature profiles as functions of 02... Here, j =1,t=0.3, Biz=0.5, Bi, =0.1 and 


6;.=0. 


Figure 4. Temperature profiles as functions of j. Here, 02..=5, 0;.=0, Biz=0.5, Bi; =0.1 and 


t=0.3. 


Figure 5. Dimensionless temperature history of the inner surface as a function of y. Here, 


21 
Biz =0.40, Bi; = 0.04, j(t)=1.14, T).=7; =300K and the T>., variation was calculated using 


FIRST [10] for a recent full scale fire [8]. 


Figure 6. Variation of the dimensionless breaking time, tT, , with dimensionless shading thickness, 


s/L for different values of absorption length, y. Here, $2 = 1.0, 61/2 = 1/3 and j/o2 =1. 


Figure 7a. Energy map representing the different heat sources and sinks for the compartment fire 


shown in Fig. 2 of Ref. 8. The parameters are the same as in Fig. 5 except here j =0. 


Figure 7b. Heat flux map representing the different heat sources and sinks for the compartment 


fire shown in Fig. 2 of Ref. 8. All parameters are the same as in Fig. 5 except here j = 0. 


Figure 1. 
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